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We determined the critical exponent v of the scalar O(N) model with a strategy based on the 
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treatment of the model shows that there is an infrared fixed point in the broken phase. The appearing 
degeneracy induces a dynamical length scale there, which can be considered as the correlation length. 
It is shown that the IR scaling behavior can account either for the Ising type phase transition in the 
3-dimensional O(N) model, or for the Kosterlitz-Thouless type scaling of the 2-dimensional 0(2) 
model. 
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I. INTRODUCTION 

The renormalization group (RG) method is a power- 
ful tool to take into account all the quantum fluctuations 
systematically [1, 2]. The RG flow equations constitute a 
bridge between the high energy, microscopic, ultra-violet 
(UV) physics and the low-energy, infrared (IR) one. One 
usually starts with a UV potential, which describes the 
small distance interactions among the elementary exci- 
tations and, by taking into account the quantum fluc- 
tuations in a consecutive manner, one finally arrives at 
the effective theory which contains the complete effects 
of the quantum modes. 

The RG flow equations are highly non-linear and, due 
to the complexity of initial conditions, it is difficult to 
map out the whole phase space, therefore one should use 
some approximations. In the vicinity of the fixed points 
of the RG equations one can linearize the flow equa- 
tions which can make the physical problem analytically 
treatable there. The 3-dimensional (3d) O(N) symmet- 
ric scalar model contains a Gaussian and a Wilson- Fisher 
(WF) fixed point. Novel approaches and new improve- 
ments of the RG method are usually tested by calculat- 
ing the scaling exponents in the vicinity of its non-trivial 
WF fixed point. The exponents can be calculated numer- 
ically by field expansion of the potential [3-5], or e.g. by 
e-expansion in 4 — e dimensions [6] , furthermore one can 
investigate the convergence the value of the exponents 
in the derivative expansion [7] . In the framework of the 
RG method more precise results can be obtained without 
expanding the potential [8, 9]. The exponents can also 
be calculated by improved Hamiltonian for the 3d Ising 
model [10]. 

Our goal is to show that there is an IR fixed point in 
the d-dimensional O(N) models inducing an IR scaling 
regime, which enables one to determine the correlation 
length and the critical exponents there. Since the IR 
fixed point is attractive, the scaling of the correlation 
length £ cannot be determined with the usual technique 
according to the linearization around the saddle-point- 
type WF fixed point. We use an non-conventional strat- 



egy in order to define £ which is based on the IR behavior 
of the RG flows in the broken phase [11, 12]. It was shown 
[13-15] that at a certain momentum scale the RG evo- 
lution stops that defines a maximal length scale which 
can characterize the appearing global condensate there. 
This scale can be identified by the correlation length £. 
It diverges as a power law-function which can provide us 
the critical exponent v of £ in the IR limit. This method 
enabled one to determine the value of v of the 3d scalar 
0(1) model, and to recover the scaling of the Kosterlitz- 
Thouless (KT) -type phase transition in the sine-Gordon 
(SG) model, too [11]. 

We use the functional RG treatment for the d-di- 
mensional O(N) model for its effective average action 
in order to get the deep IR scaling behavior of the ef- 
fective potential. We show that the exponent v based 
on the power law scaling of the correlation length in the 
IR limit for the 3d O(N) model and for the O(l) model 
with continuous dimension coincides with the one, which 
can be obtained around the WF fixed point by the con- 
ventional method based on finding negative reciprocal of 
the eigenvalue corresponding to the single relevant op- 
erator. We also show that the other critical exponents 
are not necessarily equal in these regions, e.g. the expo- 
nent rj calculated in the IR scaling regime is much larger 
than its value determined in the vicinity of the WF fixed 
point. Furthermore, we investigated the 2d 0(2) model, 
where instead of the WF fixed point a critical slowing 
down of the evolution appears. We obtained that the IR 
scaling reproduces properly the KT-type scaling of the 
model in a very simple way. It is important to note that 
our treatment can determine the value of the exponent 
v for the KT-type scaling directly, which is in principle 
not possible in the conventional way. 

The paper is organized as follows. In Sect. II we give 
the RG evolution equation of the scalar O(N) model. In 
Sect. Ill the typical phase structure of the O(N) model is 
presented and the appearing fixed points are discussed. 
The results for the exponent v for the 3d O(N) model and 
for the O(l) model with continuous dimension are shown 
in Sects. IV and V. We recover the essential scaling of the 
correlation length in the IR limit for the 2d 0(2) model 
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in Sect. VI, and finally, in Sect. VII the conclusions are 
drawn up. 



II. RENORMALIZATION 

The successive elimination of the quantum fluctuations 
is performed by means of the Wetterich RG equation for 
the effective action [1] 



1 kd k R k 

kdkTk = 2 Tr i4TTf 



(1) 



where the prime denotes the differentiation with respect 
to the field variable 0, and the trace Tr denotes the in- 
tegration over all momenta and the summation over the 
internal indices. Eq. (1) has been solved over the func- 
tional subspace defined by the ansatz 



•a + V 



(2) 



with the wavefunction renormalization Z and the poten- 
tial V, which are functions of the invariant p = cf) a 4> a /2. 
The potential has the form 



M 



2 = 2 



(3) 



with M the degree of the Taylor expansion and the di- 
mensionful couplings g n , n > 1. One can introduce the 
dimensionless couplings according to gi — k^n and 
9n = fc' 9 "'A„, with [gi] — d — 2 and [g n ] = d + n(l — d/2) 
for n > 2. For shorthand we use A2 = A. The further di- 
mensionless quantities are denoted by ~, e.g. V = k d V, 
thus 



M 



(4) 



T(d/2). In Eq. (6) we introduced the anomalous dimen- 
sion r] which is defined as r\ — —d\ogZ/d\ogk and can 
be calculated by means of the couplings as 



n 



d 1 + 2kX 



(7) 



if we consider the Litim's regulator. The inclusion of rj in 
the RG equation accounts for the evolution of the wave- 
function renormalization. We note that a more precise 
treatment can be obtained if one Taylor expand Z as the 
potential V in Eq. (4) , and considers its evolution equa- 
tion which has a similar but more involved structure as 
the potential has in Eq. (6). We also note that the in- 
clusion of the full momentum dependence of the wave- 
function renormalization [18] would make the behavior 
of the anomalous dimension regular, which could show 
that the obtained scaling of r\ is the consequence of the 
approximations used here. 

From the functional RG equation in Eq. (6) one can 
deduce evolution equations for the couplings k, A and 
A ra , n > 2, e.g. in d = 3 the flow equations are 



kdkK = —k + 
fc<9 fc A = — A + 



1 



2tt 2 (1 + 2kA) 2 

3A 2 
7t 2 (1 + 2kA) 3 



(8) 



for the first two couplings if we set 77 = and A„ = for 
n > 2. The scale k covers the momentum interval from 
the UV cutoff A to zero. In numerical calculations we 
typically set A = 1. The dependence of the number of 
couplings should always be considered [3, 4, 19]. We also 
investigated the M dependence of our numerical results, 
and we obtained that the choice M = 8 gives stable re- 
sults for the exponents in the IR limit if we Taylor expand 
the potential. 



III. THE PHASE SPACE, FIXED POINTS 



The introduction of the coupling k serves a better con- 
vergence in the broken phase. The IR regularization is 
chosen to be the Litim's regulator [16, 17] 



R k = (k 2 - q 2 )6(k 2 - q 2 ), 



(5) 



which gives fast convergence during the evolution, and 
provides simple analytic forms for the flow equations. 
The evolution equation for the potential can be derived 
from Eq. (1) [3], which reads as [16] 



kd k V 



-dv + (d-2 + n)pv' + ( 1 
1 



+ 



4«d 
d 

N - 1 



d + 2 



1 + V + 2pV" 1 + V 



(6) 



with the Litim's regulator in Eq. (5), where ' = S/6p 
and v d = l/T(d/2)2 d+1 7r d / 2 , with the Gamma function 



The O(N) model in d — 3 has two phases. The typical 
phase structure is depicted in Fig. 1 for the couplings in 
Eq. (8). There arc two fixed points in the model, which 
can be easily identified from the RG equations in Eq. (8) 
if one solves the static equations. The Gaussian fixed 
point is situated near the origin at k* g = \/2tt 2 and 
Xq = 0. The linearization of the flow equations around 
the Gaussian fixed point gives a matrix with positive 
eigenvalues (sqi — 1 and sg2 = 1) for the couplings 
showing that this fixed point is repulsive, or UV attrac- 
tive. The WF fixed point can be found at the values 
of the couplings k^f = 2/97T 2 and A^/ F = 97r 2 /8, with 
eigenvalues swfi = 1/3 and swf2 = —2. The opposite 
signs of the eigenvalues refer to the hyperbolic, or sad- 
dle point nature of the fixed point with attractive and 
repulsive directions in the phase space. 

We usually identify the critical exponent v of the cor- 
relation length £ in the vicinity of the WF fixed point 
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FIG. 1: The phase space of the 3d O(N) model. The flows 
belonging to the symmetric (broken) phase tend to right (left), 
respectively. 



by taking the negative reciprocal of the single negative 
eigenvalue, which gives t\wf = 1/2 in this case. For 
many couplings we should find the fixed points by solv- 
ing the static RG equations numerically. These equa- 
tions are non-linear, and it is very difficult to find them 
in an M dimensional phase space, especially if one does 
not know in which region of the space space they can 
be found. Conventionally the fixed points are found by 
e.g. shooting method where the non-expanded potential 
is considered. However this technique is inappropriate 
for models where only a critical slowing down appears in 
the RG flow, since no fixed points can be found. This 
situation takes place e.g. in the 2d 0(2) model. 

It is worthwhile to mention that the flows tend to a 
single curve beyond the WF fixed point. In the broken 
phase one obtains a (super)universal effective potential. 
The phase structure suggests in Fig. 1 that a further fixed 
point may exist in the IR limit. The IR effective potential 
in the broken phase can be characterized by an upside 
down parabola [13, 14] which implies that k* ir — > — oo 
and X* IR = 0, therefore the IR limit should be found 
there. Mathematically it seems that there are no fur- 
ther fixed points of the RG equation in Eq. (8), but the 
values for k* ir and X* IR could satisfy the static equations, 
although they give expressions like oo — oo and 0/0. How- 
ever if one rescales the couplings as in [20, 21] then the 
attractive IR fixed point can be uncovered. If one re- 
peats those calculation for the couplings in Eq. (8) one 
finds the following pair of evolution equations 

d T uj = 2w(l-w) - -^(3-4w) 

d T t = £(5u - 6) + —(l - u ), (9) 

where to = 1 + 2kA, i = X/uj 3 and d T — ujkdk- The static 
equations now have the Gaussian (£q — 0, oJq = 1), the 
WF (£ WF = n 2 /3, uj wf = 3/2) and the IR (£* IR = 2tt 2 /3, 



lu* ir = 0) fixed point solutions. Naturally the Gaussian 
and the WF ones has the same behavior as was obtained 
from direct calculations. However the new IR fixed point 
indeed corresponds to the values k* ir — >• — oo and A}^ = 
0, and the linearization in its vicinity gives a negative 
and a zero eigenvalue, showing that the fixed point is IR 
attractive, in accordance with the flows in Fig. 1. 

If one considers more couplings and includes the run- 
ning of the anomalous dimension one obtains qualita- 
tively similar results with similar phase space structure 
as in Fig. 1 and with the same fixed points. However, in 
the broken phase the effective potential has a wide range 
of flat region in its middle, implying that the Taylor ex- 
pansion of the potential as in Eq. (4) does not converge 
well. One can obtain more reliable results if one consid- 
ers the evolution of the non-expanded potential in Eq. (6) 
without a polynomial ansatz. Nevertheless we continue 
our investigation by expanding the potential, similarly to 
[4, 5, 16], since our goal is to demonstrate the existence 
of the IR fixed point and to calculate the corresponding 
exponents in its vicinity, and not to get better exponents 
than the ones obtained by high precision Monte-Carlo 
calculations [22, 23]. 

The RG flow for M number of couplings can also drive 
the evolution in Eqs. (6) and (7) to degeneracy where 
co = 0. The r.h.s. of the RG flow equations become sin- 
gular, when the degeneracy condition w = is satisfied at 
a certain value of k = k c . As k approaches k c the anoma- 
lous dimension r\ becomes also singular and it tends to 
zero at k c abruptly, as it is shown in the upper figure in 
Fig. 2, which also demonstrates that there are three dif- 
ferent scaling regions in its flow. In order to analyze them 
we plotted r\ as the function of the shifted scale k — k c in 
the lower figure in Fig. 2. In the UV region the Gaussian 
fixed point drives the evolution of the anomalous dimen- 
sion. In this regime it grows according to the power law 
scaling rjuv ~ {k — k c )~ 2 . There is a crossover (CO) 
scaling region between 10 ~ 8 < k — k c < 10~ 4 where a 
plateau appears giving a constant value for r\co ~ 0.043 
due to the WF fixed point. Going further in the evolu- 
tion towards the smaller values of k below k — k c ~ 10~ 8 
one can find a third scaling regime. It appears a simple 
singularity in the upper figure, but the shifted scale k — k c 
clearly uncovers the power law scaling of the anomalous 
dimension there according to m r ~ (k — kc) 1 . This scal- 
ing region is induced by the IR fixed point. 

The evolution of the other couplings also shows such 
type of scaling regimes with similar singularity structure 
in the IR limit. There the power law scaling behaviors 
also takes place as the function of k — k c with the cor- 
responding exponents. This shows that the appearing 
singularities arc not artifacts and the RG flows can be 
traced down to the value of k c . 

We note that one can find such a value of M where 
the evolution does not stop as in the upper figure in 
Fig. 2, but after the sharp fall rj continues its RG evo- 
lution marginally giving a tiny value of 77 there. How- 
ever, the singular-like fall possesses the same power-law 
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FIG. 2: The evolution of the anomalous dimension r\ is pre- 
sented as the function of the scale k (top) and the shifted 
scale k — k c (bottom). The flow of r\ has a strong singularity 
as the function of the scale k at k — k c , but it shows a power 
law like behavior as the function of the shifted scale k — k c . 



like behavior as the function of k — k c , for any value of 
M. It suggests that the value of r\ rapidly falls to zero 
at k c and it is due to the numerical inaccuracy, whether 
the RG evolution survives the falling and can be traced to 
any value of k, or the flows stop due to the appearing sin- 
gularity. It strongly suggests that the singular behavior 
with its uncovered IR scaling for the shifted scale k — k c 
is not an artifact but is of physical importance. We note 
that one has the same singularity if one solves the RG 
equation in Eq. (1) without any functional ansatz for the 
potential [9]. The RG flows stop at k c due to the huge 
amount of soft modes close to the degeneracy. It implies 
that during the evolution a dynamical momentum scale 
is generated, which can be identified by the characteristic 
scale of the global condensate in the broken phase. When 
the UV value of n\ is fine tuned to its critical value k\ 
the dynamical scale k c tends to zero. The reciprocal of 
k c can be identified by the correlation length £ of the 
model. Such an identification of £ ensures that the anal- 
ysis is performed in the vicinity of the IR fixed point. In 
the IR limit £ diverges similarly to what it does in the 
presence of any other fixed point, so the IR physics of 
the broken phase possesses an IR fixed point. 

In Fig. 3 the quantity u = 1 + 2k\ is plotted for various 
initial conditions. It shows that in the broken phase, 
where uj — > the flows break down determining the scale 
£ = l/k c as the size of the appearing condensate. In the 
symmetric phase co remains finite during the evolution 
and gives diverging flows in Fig. 3 according to power- 
law behavior in the IR. One can identify the reduced 




FIG. 3: The scaling of the degeneracy condition 1 + 2k\ as 
the function of the scale k. The diverging curves correspond 
to the symmetric phase, while the curves tending sharply to 
zero correspond to the broken phase. 
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0.990 



TABLE I: The critical exponent v in the O(N) model for 
various values of N. 

temperature t in the O(N) model as the deviation of the 
U V coupling ka to its UV critical value, i.e. t ~ k\ — k\. 
The correlation length has the power law scaling behavior 
according to 

(io) 

which characterizes a second order phase transition. The 
IR defined £ gives us a possibility to recalculate the ex- 
ponent v in the vicinity of the IR fixed point. 

IV. 3d O(N) MODEL 

High accuracy calculations exist [22, 23] for the calcu- 
lation of the exponent v in the vicinity of the WF fixed 
point for the 3d O(N) model. Our goal is now to get v 
from the scaling around the IR fixed point. We numer- 
ically calculated the scaling of the couplings and then 
we determined the critical scale k c = l/£ for given UV 
values of n\. The results are plotted in Fig. 4. For a 
given value of N we obtain power-law behavior for the 
scaling of £, and the slope of the lines provides the ex- 
ponent v in the log-log scale. The obtained results are 
listed in Table I. We denoted the WF (IR) values of v as 
I/ WF (^IR)' respectively. The results show high coinci- 
dence. The values vyfp are taken from results obtained 
from derivative expansion up to the second order, since 
this approximation are the closest to our treatment. One 
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FIG. 4: The scaling of the correlation length £ as the function FIG. 5: The value of the critical exponent v for various di- 
of the reduced temperature t, for various values of N. mension d (filled squares). The continuous curve obtained 

from e-expansion [24]. The empty squares denote the exact 

values at integer dimensions. 



can conclude that the exponent v can be also determined 
from the scaling around the IR fixed point, and has the 
same value as was obtained around the WF fixed point. 
This fact might seem quite strange at a first glance since 
the two fixed points are well separated in the phase space 
according to Fig. 1. The coincidence may come from the 
fact that the condensate of the broken phase has a global 
feature which accompanies the whole flow starting from 
UV to IR. The other exponents do not necessarily coin- 
cide, e.g. anomalous dimension rj « 0.04 in the vicinity 
of the WF fixed point, while it converges to zero in the 
IR one. 



V. DIMENSION-DEPENDENCE OF v IN THE 
O(l) MODEL 



gle in Fig. 5, which does not appear in [26]. This is due 
to the numerical inaccuracy in our method coming from 
the high anomalous dimension in the IR regime which 
should be traced to extremely small values of the scale 
k. We expect that by taking into account higher orders 
in the derivative expansion or by treating the problem 
with full momentum dependence [18] the wiggles could 
be removed. 

These results clearly show that it is possible to deter- 
mine the power-law scaling of the correlation length in 
the vicinity of the IR fixed point. We note that there 
is no WF fixed point in case of d — 4 model. However 
the degeneracy induced scaling gives us back the analytic 
result even in this case, too. 



As a further test we calculated the exponent v for 
the O(l) model for continuous dimension. By using e- 
cxpansion the dimension dependence of v can be easily 
obtained [24, 25]. By linearizing the functional RG flow 
equations around the WF fixed point the exponent v can 
be calculated for any dimensions [26] . Naturally far from 
d = 4 the anomalous dimension grows up, which requires 
to take into account the evolution of the wavefunction 
renormalization. 

We determined the exponent v by the degeneracy in- 
duced scaling around the IR fixed point. The results are 
shown in Fig. 5. The exponent v calculated around either 
the WF or the IR fixed points give the same results. The 
square at d = 2 in Fig. 5 shows v = 1, as the exact results 
of the 2d Ising model [27] . When d — 3 there are no exact 
calculations but one can find a huge amount of articles 
related to calculating v [22, 23]. When d = 4 the mean 
field calculations give v — 0.5. Our results give better re- 
sults than the e-expansion ones for low dimensions, since 
far from d = 4 one needs more and more loop corrections 
to get reliable results. At around d « 2.5 there is a wig- 



VI. 2D 0(2) MODEL 

The 2-dimensional 0(2) model is exceptional among 
the O(N) models since it possesses an infinite order KT- 
type phase transition [28]. This transition is character- 
ized by the scaling of the correlation length according 
to 

log^t"", (11) 

with the exponent v — 0.5. The 2d 0(2) model belongs 
to the same universality class as the 2d XY model [27], 
the 2d Coulomb gas [29] the 2d SG model [30], but the 
equivalent models inherently contain vortices as elemen- 
tary excitations which can account for the essential scal- 
ing. The functional RG approach shows KT-type scaling 
according to Eq. (11) if one considers the effect of the 
wavefunction renormalization in the SG model [21], fur- 
thermore it was found that the KT and the IR scaling 
gives the same scaling behavior and the same exponent 
v « 0.5 [11]. However the 2d 0(2) model does not contain 
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FIG. 6: The KT-type scaling of the correlation length in the 
0(2) model in d — 2. The lines correspond to flows belonging 
to different UV values of A. 



any vortices but the scaling in Eq. (11) can be uncovered 
[31, 32], although in a rather involved way. 

The 2d 0(2) model also contains an IR fixed point, 
therefore the degeneracy induced IR scaling also can be 
used. We plotted the results in Fig. 6. They show that 
we can recover the essential scaling according to Eq. (11) 
in a very simple way by using our proposed technique, 
although the exponent shows a slight dependence on the 
UV value of A as is depicted in the inset of Fig. 6. Besides, 
one can conclude that v « 0.5, in accordance with other 
results [21, 31, 33]. 

VII. SUMMARY 

The critical exponent v of the correlation length £ 
was calculated for the O(N) model in the IR limit. We 
showed that there is an IR fixed point in the broken phase 
of the model, which induces a degeneracy and stops the 
RG evolution at a finite momentum scale. The latter de- 
fines a characteristic length scale which can be identified 
with the correlation length. This technique of identify- 
ing the correlation length enabled us to determine the 
exponent v around the IR fixed point, far from the WF 
one in which vicinity it is usually calculated. We showed 
that the value of the exponent v agrees well with results 
obtained in the vicinity of the WF fixed point for the 



3d O(N) model, see Table I. The qualitative behavior 
of v for the 0(1) model with continuous dimension pre- 
sented earlier by e-expansion or linearization around the 
WF fixed point was recovered by the IR scaling. Further- 
more the essential scaling of the 2d 0(2) model with the 
exponent v ss 1/2 was also got back. The latter model 
exhibits a KT-type phase transition showing that this 
simple technique can be easily applied to describe any 
types of phase transitions. 

The coincidence of the exponent v around the WF and 
the IR fixed point may come from the fact that the corre- 
lation length in the broken phase characterizes the con- 
densate which occurs there and consists of a macroscopic 
population of soft modes showing the global nature of 
£. The anomalous dimension r\ does not agree around 
the WF and the IR fixed points. The common feature of 
the d-dimensional O(N) models for d > 2 that there is a 
crossover fixed point, namely the WF fixed point there. 
Although the crossover and the IR fixed points are far 
from each other in the phase space, we read off the scal- 
ing of £ from the same trajectories which first approach 
the crossover fixed point and then tend towards the IR 
one. However the real power of this technique can emerge 
if there are no crossover scalings. For example in the 2d 
0(2) model there is only some remainder of the WF fixed 
point as a critical slowing down of the RG flows, or in 
the 4d O(l) model the WF fixed point melts into to the 
Gaussian one. It is shown that the IR scaling can cor- 
rectly determine the order of the phase transitions, and 
the proper value of the exponents v in these models, too. 
A further non-trivial example is the massive SG model 
[34], and the layered SG model [11, 35] where there are 
no crossover fixed points too, but the IR scaling uncovers 
the type of the phase transition and helps us to determine 
in which universality class the models belong to [36] . 
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